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SUMMARY 
New  error  bounds  are  developed  for  the  Liouville- 
Green  approximation  to  the  solution  of  an  important  class 
of  differential  equations  arising  in  military  operations 
research  (specifically,  variable-coefficient  Lanchester- 
type  equations  of  modern  warfare  for  combat  between  two 
homogeneous  forces) .   In  contrast  to  many  previous  re- 
sults, our  error  bounds  apply  to  initial-value  problems 
and  are  expressed  in  terms  of  initial  conditions.   Previous 
error  bounds  for  boundary-value  problems  are  sharpened  as 
a  consequence  of  our  development  of  these  new  error  bounds 
for  initial-value  problems.   Finally,  applications  are 
made  to  some  important  specific  models  of  combat  between 
two  homogeneous  forces  with  time-dependent  attrition- 
rate  coefficients. 


1.   Introduction 

OLVER  [12]  has  developed  error  bounds  for  the  so-called 
LIOUVILLE-GREEN  '     (LG)  approximation  [5,  10]  to  the  solution  of 
the  differential  equation 

d2x 


dt 


7    =  J(t)x.  (1) 


His  results  extended  earlier  work  by  BLUMENTHAL  [1] .   The  LG 
approximation  is  of  particular  importance  in  applied  mathematics 
because  it  involves  only  elementary  functions.   As  OLVER  [12] 
has  pointed  out,  however,  the  development  of  strict  upper  bounds 
for  the  errors  in  the  approximate  solutions  has  received  rela- 
tively little  attention.   Moreover,  the  error  bounds  that  have 
been  previously  developed  [12,  13]  are  for  boundary-value  prob- 
lems, since  most  problems  of  interest  in  mathematical  physics 
are  boundary-value  problems .   These  error  bounds  do  not  apply 
to  initial-value  problems. 

Thus,  the  purpose  of  this  paper  is  to  develop  error  bounds 
for  the  LG  approximation  to  initial-value  problems.   Furthermore, 
in  recasting  OLVER' s  results  in  a  form  suitable  for  initial-value 
problems,  we  have  been  able  to  sharpen  his  bounds  for  boundary- 
value  problems.   Although  we  will  develop  our  results  within 
the  context  of  a  specific  problem  in  operations  research,  they 
are  clearly  applicable  to  the  general  second  order  initial-value 
problem,  and  it  is  a  straightforward  task  to  so  recast  them. 


Also  called  the  WKB  approximation  [13]. 


Probably  the  most  widely  used  (at  least  in  the  United 

States)  deterministic  differential-equation  model  in  operations 

2) 
research  [4,15-17]  are  LANCHESTER-type  equations  of  warfare   , 

which  are  named  for  the  pioneering  1914  work  of  F.  W.  LANCHESTER 
[9],   In  this  paper  we  consider  such  a  linear,  variable-coeffi- 
cient differential-equation  model  for  combat  between  two  homogene- 
ous forces.   This  model  yields  an   X   force-level  equation  equiva- 
lent to  (1).   Unfortunately,  for  even  the  simplest  time-varying 
attrition-rate  coefficients  of  interest,  the  solutions  cannot 
in  general  be  expressed  in  terms  of  either  elementary  functions 
or  tabulated  higher  transcendental  functions  [17] .   It  is  there- 
fore natural  to  seek  an  approximation  in  terms  of  elementary 
functions.   Moreover,  one  is  always  interested  in  a  simple 
a  priori  estimate  for  the  error  in  the  approximate  solution. 
In  this  paper  we  develop  an  error  bound  that  is  both  realistic 
and  easy  to  evaluate.   This  bound  is  then  computed  for  some 
particular  attrition-rate  coefficients  of  interest. 

2.   Variable-Coefficient  LANCHESTER-Type  Equations  of  Modern 
Warfare 

We  consider  the  following  variable-coefficient  LANCHESTER- 
type  equations  of  modern  warfare    for  combat  between  two 
homogeneous  forces 


2) 

We  will  refer  to  any  differential-equation  model  of  combat  as 

being  LANCHESTER-type  equations.   The  state  variables  are 

typically  the  force  levels  of  the  various  weapon  system  types. 

3) 

The  term  "of  modern  warfare"  denotes  that  we  are  considering 

linear  differential  equations.   There  are  other  nonlinear  types 

of  LANCHESTER  equations  [4,16]. 


dx/dt  =  -a(t)y/  dy/dt  =  -b(t)x,        (2) 

with  initial  conditions 

x(t=0)  =  xQ,  y(t=0)  =  yQ 

where   t=0   denotes  the  time  at  which  the  battle  begins,   x(t) 
and  y(t)   denote  the  numbers  of   X   and   Y   at  time   t,   and 
a(t)   and  b(t)   denote  time-dependent  LANCHESTER  attrition- 
rate  coefficients.   These  equations  (2)  have  been  hypothesized 
to  model  combat  in  which  both  sides  use  aimed  fire  and  target 
acquisition  times  are  independent  of  the  numbers  of  firers  and 
targets  [2,19].   The  attrition-rate  coefficients  represent  the 
effectiveness  of  each  side's  fire  (i.e.  its  firepower).   Temporal 
variations  in  a  side's  fire  effectiveness  (caused  by,  for  example, 
changes  in  force  separation,  combatant  postures,  target  acquisi- 
tion rates,  etc.)  are  modelled  by  the  time-dependent  attrition- 
rate  coefficients.   Further  discussions  of  the  physical  assump- 
tions hypothesized  to  yield  (2) ,  estimation  of  the  attrition- 
rate  coefficients,  and  the  importance  of  (2)  to  military  opera- 
tions research  are  found  in  [16,17]. 

We  assume  that   a(t)   and  b(t)   are  positive  and  twice 
differentiable  for   tfl  <  t  <  +°°  with   tQ  <  0 .   We  also  assume 
that   a(t),  b(t)  £L(tQ,T)   for  any  finite   T.   We  further  take 

a(t)   and   b(t)   to  be  given  in  the  form   a(t)  =  k  g(t), 

a 

b(t)  =  k,h(t),   where   k  ,  k,   are  positive  constants  chosen  so 
b  a   b 

that   a(t)/b(t)  =  k  /k,   when   g(t)  =  h(t)   for  all   t.   We 

a   d 

introduce  (see  [17])  the  intensity  of  combat   I(t)   and  the 
relative  fire  effectiveness   R(t)   defined  by 


I(t)  =  /a(t)b(t)  ,      and     R(t)  =  a(t)/b(t).      (3) 


We 


similarly  introduce  the  combat-intensity  parameter   X_   and 


the  relative- fire-effectiveness  parameter   A    defined  by 

"" — — — — — — — — — — — — — — — — — — — — — — — — — — — — —     i\ 

X_  =  /kTkT,        and        XD  =  k  /k.  .        (4) 
lab  R    a  b 

A  large  class  of  tactical  situations  of  interest  can  be 
modelled  with  the  following  general  power  attrition-rate 
coefficients  [17] 

a(t)  =  k  (t+C)y,       and    b(t)  =  k,  (t+C+A)V/      (5) 
a  d 

where   A,C  >   0.   We  will  call   A   the  offset  parameter,  since  it 

allows  us  to  model  (with   y,v>0)   battles  between  weapon  systems 

with  different  maximum  effective  ranges.   We  will  call   C   the 

starting  parameter ,  since  it  allows  us  to  model  (again,  with 

u,v>0)  battles  that  begin  within  the  minimum  of  the  maximum 

effective  ranges  of  the  two  systems.   The  offset  and  starting 

parameters  are  related  to  various  physical  quantities  in  [17] . 

We  observe  that   t  =-C.   Also,   a(t),  b(t)  £L(t  ,T)   implies 

that   y,v  >  -1. 

From  (2)  we  may  obtain  the   X   force-level  equation 


d2x 


-{irErat}at-a<t)b<t>x  =  0'  (6) 


dt 
with  initial  conditions 

x(t=0)  =  xQ,        and   { [1/a (t) ] dx/dt }t=Q  =  -yQ. 

X   Y  X 

We  may  consider  that   tQ  =  max(tQ/t0)/   where   tQ   denotes  the 

right-most  finite  singularity  of  the   X   force-level  equation. 

Furthermore,  we  set   ^  =  0   if  there  are  no  finite  singularities 


3.   LIOUVILLE-GREEN  Approximation  to  LANCHESTER-Type  Equations 
of  Modern  Warfare 

Let 


t  =      /a(s)b(s)  ds  = 
fc0 


I(s)ds,  (7) 


and  denote   x(t=0)   as   x  .   Then   x  >  0   for   t  <0.   From 
a(t),  b(t)  £L(t0/T)   it  follows  that   x  =  x(t)   is  well  defined 
by  the  CAUCHY-SCHWARZ  inequality  for  integrals.   The  transforma- 
tion is  also  easily  seen  to  be  invertable .   We  observe  that 
t  -  xn   is  related  to  the  average  intensity  of  combat   I(t)   by 


*  ||  J   I(s)ds|t  =  tl(t)  . 


We  may  call  T  -  x    "the  elapsed  normalized  battle  time/"  since 
the  transformation  (7)  reparameterizes  the  battle's  time  scale 
in  terms  of  elapsed  time  of  battle  and  average  combat  intensity 
as  (8)  shows  us.   The  substitution  (7)  transforms  (6)  into 

d2x   jl  d  „  „/.  x)  dx       -  ,nv 

with   initial    conditions 

x(x=xQ)    =   xQ/      and      { [1/R1/2 (t) ] dx/dx }t=t      =   -yQ . 

Remark    1:         It   is   easily   shown    (see    [17])    that    (9)    may  be 
transformed  into  a   linear  second  order  differential  equation 
with   constant   coefficients   if  and  only   if 

1        d 


I(t)    dt 


In  R(t)    =   constant. 


Let   Rn   denote   R(t=0) .   The  substitution 

x(t)  =  X(t) [R(t)/RQ]1/4,  (10) 

transforms  (9)  into  LIOUVILLE's  normal  form  (see  INCE  [7]  or 
KAMKE  [ 8 ] ) 


d2X 


dT 
with  initial  conditions 


?-  U+F(t)  }X  =  0,  (11) 


X(x=xQ)  =  xQ,   and   dX/dx(x=xQ)  =  -yQfR^-  xQ£Qt 
where 

F(t)  =  P"(t)/P(t) ,  P(x)  =  [R(t)]~1/4,      (12) 


(13) 


£(t)  =  JirtTdT*11*' 

e    denotes   e (t=0)  ,   and   P*(t)   denotes   dP/dx . 

Writing  (11)  as   d2X/dx2  -  X  =  F(x)X,   we  may  use  variation 
of  parameters  to  obtain  the  solution  to  (11)  as 

X(t)  =  xQ  cosh  (t-Tq)  -  (y0/R^+x0e0)sinh(T-T0) 

-X 

+    F(a)sinh  (x-a)X(a)da.     (14) 


F(a; 


0 

4) 
If  one  drops    the  integral  term  in  (14) ,    one  obtains  the 

LIOUVILLE-GREEN  approximation   X(x) 


4) 

Heuristically,  if  the  appropriate  fractional  power  of  the  relati 

effectiveness   R(t)   is  "slowly  varying,"  then  by  (12)  we  would 

expect  that   |F(x) |<<1   so  that  the  integral  term  in  (14)  is 

"negligible."   Theorem  1  gives  us  bounds  on  how  "negligible" 

this  term  is. 


X(t)  =  xQ  cosh  (t-t0)  -  (y0y/R^+x()e0)sinh(T-T0)  ,       (15) 
which  in  terms  of  the  original  independent  variable   x   reads 


x(t)  = 


*&\         {xQcosh  (t-t0)  -  <y0/R^-xo£o)sinh(T-T0)}.  (16) 


We  observe  from  (14)  that   F(t)  >0   for  all   t  >  T-.   implies  that 

as  long  as   x(t)  >0  we  have   x(t)  >x(t)  .   A  similar  statement 

holds  for   F(t)  5  0.  As  we  shall  see  below,  such  cases  in  which 

F(t)   is  always  >0  or  5  0   are  readily  encountered  in  applica- 
tions . 

4 .   Error  Bounds  for  the  LIOUVILLE-GREEN  Approximation 
The  main  result  of  this  paper  is  Theorem  1. 

Theorem  1:   Error  bounds  for  the  LIOUVILLE-GREEN  approximation 
are  given  by 

|x(t)-x(t)  |  <x.  K_e(t)    <  x.  K     e(t)  ,  (17) 

J     j  j     u 

where 

y 

Ky  =   2t(l+|e0l)  +^^>  ,  (18) 

j  =  I      for      1-|£^<e0       and     Kj  =  l+e0  +  ^^,  (19) 

y  y 

j  -U      for      -l-_i^<£()<l--J^       and      K„   -   2.         (20) 

J  =   III      for      ^-1-^^      and     K^   =  1  -  eQ  -  ^  ^>  0,       (21) 
and 


e(t)  = 


R(t) 


nl/4 


«-R0  - 


{exp  (• 


ri 


iQ 


F(a) |da)-l}sinh(x-T  ) . 

0 


(22) 


The  sign  of  the  error  is  determined  by  the  sign  of   F(t) .   As 
long  as   x(t)  >0,   it  follows  that 

F(t)  >  0   for  all   t  >  0   implies  that   x(t)  >  x(t)  , 

with  the  last  inequality  being  reversed  when   F(t)  5  0. 

Proof:   Theorem  1  readily  follows  from  Lemma  2,  which  is  proven  below. 

5 .   Development  of  Error  Bounds 

Consider  the  following  fundamental  system  of  solutions 
{x1/x2>   to  (11) 

Xk(i)  =  {l+hk(T)}exp[(~l)(k"1)  (t-tq)  ]   for   k  =  1,2,  (23) 

where   ni(T)   for  k  =  1,2   is  to  be  chosen  so  that 


hk(T=T0)   =  0, 


and     dh  /dT(x=TQ)  =  0. 


(24) 


It  follows  that  the  solution  to  (11)  may  be  expressed  as 


1  (T_Tn) 

X(x)  =  ^-{(x0(l-e0)-y0^)e     ° 


[l+h1(x)J 


-<T-T  ) 

+  (x0(l+e0)+y0/R^)e     u  [l+h2(T)]>, 


(25 


so  that 


X(t)-X(t)  =  ^{[x0(l-£0)-y0/R^]e     u  h1(T) 


+    [x0(l+e0)+y0^]e 


-(T-T     ) 


h2(T)}. 


(26) 


bstituting  (23)  into  (11),  we  find  that   ^,(1)   for   k  =  1,2 
st  satisfy 

dX       k   dhk 

-%?£-   (-DK  2-^-F(T)hk  =  F(T),  (27) 

th  initial  conditions  (24) . 

We  may  consider   ni,(T)   f°r   k  =  1,2   to  be  an  error  term 
r  the  LI OUVILLE- GREEN  approximation.   We  next  develop  a  bound 
its  magnitude,  which  sharpens  earlier  results  by  OLVER  [12]. 

mma  1 :   A  bound  on  the  magnitude  of   nT-(T)   for   k  =  1,2   is 
ven  by 


t)  |exp  [(-l)(k"1)(x-T  )  ]  <  2(exp(i 


0 


I F (a) |da)-l}sinh(x-T0) .   (28) 
T0 


oof:      For  notational   convenience,   we   develop   the  bounds    for     h    (t)      and 
(t)      separately.      Transposing   the   right-most   terra  on   the   left-hand  side   of 
7)    for     k  =   1,      treating      {l+h    (t)}f(t)      as   a   "forcing  term,"   and   inte- 
iting  twice;   we   obtain   the   following  VOLTERRA  integral  equation   after  a 
rther   integration  by  parts 


VT,     -    i 


/•T 

{l-e2(0"T)}F(a){l+h1(a)  }do.  (29) 

To 


Lving  (29)  in  the  usual  manner  by  successive  approximations,  we  obtain 

00 

h.  (t)  =  I     T  (t),  (30) 

1       L.      n 
n=l 

ire     T  (x)  =  1  and  for  n>l 


T  C,)  -  i 


'T       _  ,     . 

{l-e^ta"T;}F(a)T  .(a)da.  (31) 

n-1 

T0 


2(a-T)        "2(T_T0) 
Observing  that   {l-e      }  <  1  -  e  for  0  <  t  -  a  <  t  ,   we  find  that 


j         -2(t-t  )\ 

l^tx)  I  5  ^"e  2 '-       |F(o)|do, 


with  equality  holding  for   t  =  T  .   A  straightforward  inductive  argument 
along  the  usual  lines  now  yields 


La-2'-*'} 


|T  (T)|5liZ£_^ L;j   |P(a)|do|   .  (32) 

2   n!      (Jx0        ) 

One   step   in  the   inductive  proof  of    (32)    is   deserving  of   further  elaboration, 

however.      From    (31)    and    (32)    we   obtain 

i  fT  oi         \     i        "2(o-T    )^  ,  r<3  x  n 

lTn+l(T)l-"nTl U-e     a"T'}|l-e  °}|f(o>|{I       |p(u)|du|do.         (33 

2  ni      tq  TQ 

The   inductive  proof  of    (32)    is   completed  by   combining    (33)    with   the  observa- 


tion  that      {l-e  )^l-e  °/<l-e  0<t„5o5t.      Our 


-2(o--xn)\  _2(t-tJ 

for      0  <  t. 
0 

sharpening  of  OLVER's    results   is   due   to   this   observation.      The  remaining 

steps   in   the  proof  of    (28)    for     k  =   1      follow  along  well-known  lines    [12,13] 

and  will  be  omitted  here.      Similar  arguments   are   used  to  prove  (28)    for 
k  =   2. 

Remark    2  :       In  our  notation,   OLVER's    [12]    corresponding 
error  bound   for    k  =  1    would  read 

|h    (t)  I    ^exp(|-J         |F(a)|do)-l    .  (34) 

2      T0 

His  corresponding  error  bound  for  k  =  2  is  not  directly 
comparable  to  our  result  here,  since  he  does  not  take 
both  errors  zero  at  the  same  point. 


5) 


In  our  notation  for  the  finite  or  infinite  interval   ^Tn,Ti^' 
OLVER  [12,13]  takes   h  (t=tq)  =  dh^di  (t  =t  Q)    =  0   and 
h2(T=T  )  =  dh  /dT(x=x  )  =0   in  contrast  to  (24).   This  is 
what  makes  his  results  unsuitable  for  initial-value  problems. 

l  n 


Remark  3:  Similar  arguments  may  be  used  to  develop  a 
bound  on  |h'(x)|.  For  our  applications  this  result  is 
not   important. 

Lemma    2 :       In    terms    of    the    transformed   dependent    variable      X(x)/ 
error   bounds    for    the   LIQUVILLE-GREEN    approximation    are    given   by 

|X(t)-X(t)|<x()KjE(t)   <xqKuE(t)/  (35) 


whe 


re      K        and      K        are    given   by    (18)    through    (21) ,    and 


E(t)    =    (exp(| 


|F(a) |do)-l}sinh(T-T0) .  (36) 

T0 


Proof:    From  (26)  and  Lemma  1  we  obtain 

|x(t)-x(t)  |  <  {|x0(l-e0)-y0^"|  +  |x0U+ea)+y0*^|}E(T).         (37) 
It  follows  that  a  rather  loose  error  bound  is  given  by 

|x(T)-X(T)  |  <xQK    E(t)  .  (38) 

^0     / —  ^0     i — 

We  observe   that   for     -1 /R„  <  e    <  1 *R„,      we   have 

xQ        0        0  xQ        0 

lxo(1-£o)-yo^l  +  lxo(1+eo)+yo/^'l  =  2V  (39) 

so   that    (37)    becomes      |  X(t)  -X(t)  |  <  2x  E  (t)  .      Thus,    (35)    is  proven   for     J  -   II 

For     1-    —  i/rT<  e    ,      the  error  bound    (37)    becomes 
x0  ° 

yn    , 

|x(t)-X(t)|  <2x.(e.+-ii  ^T)E(t),  (40) 

'00x0 


since 


VW^tAI  +  IV1+£o)+y<AI  =  2Veo+^  S*j-  <41) 
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The   error  bound    (40)    may  be   sharpened,    however,    as    follows.      If     F(a)  >0      fo] 
0<t    <a5x,      then   as    long  as      X(t)>0     we  have      X{t)>X(t)      and     g    (t)>0, 


where      for     k  =   1,2 


gk(T)    =  Jhk(^)  exp[(-l)k(T-T0)].  (42) 


y 

It   follows   from    (26)    that  for     1-—  /r~<  en 

xQ        0        0 


0<X(t)  -  X(x)  <  (x0(l+e0)  +yQ^>E+(T)  ,  (43) 


where 

ft 


E    (x)    =    {exp(r- 


F(a)da)  -  l}sinh(T-T   ).  (44) 

T0 


Similarly,    if     F(o)$0      for     t    SoSt,      then  as   long  as      X(t)>0     we  have 
X(t)<X(t)      and     g    (t)S0,      so  that 

0>  X(x)  -  X(t)  >  {x0(l+e0)+y0^}E_(x)  /  <45> 


where 


E_(t)    =   {l-exp(-j 


It   follows   from    (43)    and    (45)    that 


F(a)da)}sinh(T-TQ) .  (46) 

T0 


|x(t)-X(t)|    <  x„{l  +  e    +  —  /r~}E(t).  (47) 

1  '00x0 

y  y  y 

Furthermore,    since     e    >  1 SrI     implies   that     2(e    + —  /r~)  >  l  +  e    + —  /r~, 

0x0*  0x0  0x0 

the  bound  given  by    (47)    is   sharper  than   that  given  by    (40).      Thus,    (35)    is 

proven   for     J  =   I.      The  proof  of    (35)    for     J  =  III   is   similar   to   that   for 

J  =   II.      It  is  easily  seen   that     K    <K        for     x^,y    >  0. 

J        U  0^0 
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6 .   Examples 

We  now  compute  theoretical  error  bounds  for  two  special 
cases  of  the  LIOUVILLE-GREEN  approximation  to  the  solution  of 
(6)  with  the  general  power  attrition-rate  coefficients  (5): 
(I)  power  attrition-rate  coefficients  with  no  offset  (i.e.   A  =  Q), 
and  (II)  linear  attrition-rate  coefficients  with  positive  offset 
(i.e.   A  >  0)  . 

6.1.   Power  Attrition-Rate  Coefficients  with  No  Offset 

In  this  case  we  have 

a(t)  =  ka(t+C)u/       and      b(t)  =  kb(t+C)v,      (48) 

with   C>0   and   p,v > -1.   The  LIOUVILLE-GREEN  approximation  to 
the   X   force  level  is 

kit)    =    (1+t/C)  (y"v)/4{x0cosh  (t-t0) 

-[y0/X^  c(u-v)/2  +  Xq4^~V   c"6]sinh(T-T0)},   (49) 


where 


and 


T(t)  =  (l/6)XI(t+C)6/  (50) 


6  =  (y+v+2)/2.  (51) 


In  preparation  for  estimating  the  error  in  the  LIOUVILLE- 
GREEN  approximation  (49)  by  Theorem  1,  we  compute 

F(T)  =  (u-v)(3^v+4)  m  (52) 

166^x^ 
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Thus,   F(t)  >  0   for  all   t  >  x  >  0   if  and  only  if   y >  v.   Hence, 
as  noted  at  the  end  of  Section  3,  one  frequently  encounters  in 
applications  cases  in  which   F(t)   always  has  the  same  sign. 
For  the  error  estimate  (17) ,  we  then  have 


1 
2 


|F(o,  |d0    =   I^i|3y±v±4I{c-5_(t+cr6K  (53) 

T0  Z 


Remark  4:    The  exact  solution   x(t)   to  (6)  with  attrition- 
rate  coefficients  (48)  is  given  in  TAYLOR  and  BROWN  [17]  (see 
also  [16]).   It  may  be  expressed  in  terms  of  modified  Bessel 
functions  of  the  first  kind  of  (for   p,v>-l)  fractional  order, 

i.e.   I    for  0  <  ct  <  1.   Since  few  of  the  latter  are  tabulated 
a 

(i.e.  tabulations  only  exist  for  a  =  ±  1/4,  ±  1/3,  ±  1/2,  ±  2/3, 
±3/4,   and  these  do  correspond  to  cases  of  interest),  TAYLOR 

and  BROWN  [17]  suggested  the  use  of  new  transcendents  which 

6) 
they  called  LANCHESTER-CLIFFORD-SCHLAFLI  functions   . 

6.2.   Linear  Attrition-Rate  Coefficients  with  Positive 
Offset 

a(t)  =  k  (t+C)  ,       and      b(t)  =  k,  (t+C+A),     (54) 
a  b 

with   A,C>0.   The  LIOUVILLE-GREEN  approximation  to  the   X   force 
level  is 


After  earlier  work  by  W.  K.  CLIFFORD  [3]  and  L.  SCHLAFLI  [14] 
(see  also  [6,18]  )  . 
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~tJm*     r  (i+A/c)  i1/4r       u  ,      , 

x(t)  =  LTiWItTcTTj       {x0cosh  (t-t0) 


yo/x^+  xoA/c 


./I+A/C      4XTC2(1+A/C)  3/2_ 


sinh(x-T0)},         (55) 


where 


and 


A2  ^—  r-r- 

x(t)    -  ^-  Ajdfufip^T-  ln(i[»+/ip=T)}#  (56) 


ip  (t)    =    1  +  2(t+C)/A.  (57) 


For  estimating  the  error  in  the  LIOUVILLE-GREEN  approxi- 
mation (55),  it  is  more  convenient  to  express   F(t)   as  defined 

by  (12)  in  terms  of  the  original  independent  variable   t.   We 

7) 
compute  that 

,(T)  =    A{12(t+C)+7A}    >  (58) 

16Xj(t+C)  (t+C+A) 

It  follows  that   F(t)  >  0   for  all   x  >  x0  >  0.   Thus,  as  in  the 
previous  example,   F(t)   always  has  the  same  sign.   For  the 
error  estimate  (17)  ,  we  then  have 

0  <  j     F(a)da  <  min(m1(t) ,m2(t) ) ,  (59) 

T0 


where 


m,(t)  =  - 3A  ,/?         (4  (l-q1/2)  +  ^(l-q3/2)  }  ,     (60) 

1  16C  (1+A/C)  '  Aj               4C 

m  (t)  =  V-  (4(l-q3)  +]^(l-q4)},  (61) 

z  32CJXI            4C 


7>T  1 

In  general 


F(T>  =^r^{-^lnb(t,-7^4nR(t)+^'ntt} 
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and 

q(t)    =    l/(l+t/C)  .  (62) 

It  may   be   shown   that      m    (t,)    =   m    (t,)       implies    that      m,  (t)   >m    (t) 
for   all      t  >  t,  .      The   error   control    term  estimate      m..  (t)      was 
developed    for    "small"      t8)  ,      while      m2 (t)       for    "large"      t9) . 

Remark    5 :        This   case   is   of  particular   interest   in  military 
operations   research,    since   it  may  be   used  to  study   combat 
between   two  weapon   systems  with   different  maximum  effective 
ranges    [16].      The   exact   solution     x(t)      to    (6)      with   attrition- 
rate   coefficients    (54)    is   given   in    [17]    (see   also    [16]).      It 
apparently  cannot  be   expressed  in   terms  of  previously   "known" 
transcendents,    since   the      X      force- level   equation   in   this   case 
could  not  be    found  to  correspond  to  any  second  order   linear 
equation   considered   in    [8]    or    [11]. 

7 .      Final    Remarks 

Although  given  within  the  context  of  a  specific  problem 
in  operations  research,  the  reader  will  have  no  trouble  trans- 
lating  the   above   results    into   those    for   the   general    second  order 


8).^.             T          f1         {12(s+C)+7A}ds  -.     -    ..      . 

Letting      J   =  b f-t-= — ' - — cyy,      we    fmd   that 

J0  (s+C)*/Z(s+C+A)*/Z 


Here  we  have  used  the  bound   J  < 


F(a)do  = 
^{12  (s+C)+7A} 


T0 


16XI  (C+A)5 

9) 


0    (s+C)5/" 


In  this  case  (to  be  contrasted  with  the  previous  one) ,  we  have 

t 


used   the    bound      J  < 


{12 (s+C)+7A}ds 
(s+C)5 
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initial-value  problem.   In  the  examples  of  Section  6  we  saw 
that  for  two  models  of  considerable  importance  in  military  opera- 
tions research  both  the  LIOUVILLE-GREEN  approximation  and  bounds 
on  its  error  were  simply  expressed  in  terms  of  elementary  func- 
tions.  No  previous  application  of  the  LIOUVILLE-GREEN  approxi- 
mation has  appeared  in  the  operations  research  literature.   In 
a  subsequent  paper  we  plan  to  present  a  numerical  investigation 
of  the  accuracy  (both  numerical  eveluation  of  the  theoretical 
error  bounds  and  a  comparison  with  the  exact  solution  when 
available)  of  the  LIOUVILLE-GREEN  approximation  to  the  solution 
of  LANCHESTER-type  equations  of  modern  warfare  for  combat  between 
two  homogeneous  forces. 
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